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f -s  Abstract 

A  composite  grid  was  generated  In  an  attempt 
to  Improve  grid  quality  For  a  typical  turbine  blade 
with  large  camber  In  terms  of  mesh  control,  smooth¬ 
ness,  and  orthogonal Hy .  This  lumpostte'^grld  con¬ 
sists  of  the  C-grld  -fiwMl-grttffc'  In  the  Immediate 
vicinity  of  the  blade  and  the  H-grld  In  the 
upstream  region  and  In  the  middle  of  the  blade 
passage  between  the  C-grlds.  It  provides  ar^good 
boundary  layer  resolution  around  the  leading-edge 
region  for  viscous  calculation,  has  orthogonality 
at  the  blade  surface  and  slope  continuity  at  the 
C-H  (or  O-H^?  Interface,  and  has  flexibility  In 
control  ling  4S^mesh  distribution  In  the  upstream 
region  without  using  excessive  grid  points.  This 
composite  grid  eliminates  the  undesirable  qualities 
of  a  single  grid  when  generated  for  a  typical  tur¬ 
bine  geometry. 

A  finite-volume  lower-upper  (LU)  Implicit 
scheme  can  be  used  In  solving  for  the  turbine  flows 
on  the  composite  grid..  This  grid  has  a  special 
grid  node  that  Is  connected  to  more  than  four 
neighboring  nodes  In  two  dimensions  and  to  more 
than  six  nodes  In  three  dimensions.  But  the 
finite-volume  approach  poses  no  problem  at  the 
special  point  because  each  Interior  cell  has  only 
four  neighboring  cells  In  two  dimensions  and  only 
six  cells  In  three  dimensions.  The  finite-volume 
Ul  -Implicit  , scheme  was  demonstrated  to  be  robust 
and  efficient  for  both  external  and  Internal  flows 
In  a  broad  flow  regime.  *- 

Introduction 

Turbine  blades  are  often  designed  to  have  sub¬ 
stantial  thickness  and  camber  and  rounded  leading 
and  trailing  edges.  For  a  typical  turbine  blade 
the  H-grld  does  not  provide  a  good  boundary  layer 
resolution  around  the  leading  edge.  If  very  fine 
meshes  are  used  to  resolve  this,  the  number  of  grid 
points  In  the  upstream  region  becomes  excessive. 
Thus  mesh  points  would  be  wasted.  A  turbine  flow 
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analysis  on  an  H-grld  Indicated  that  the  standard 
H-grld  needs  cusps  at  the  leading  and  trailing 
edges  to  reach  a  convergent  solution  and  that  It 
causes  somewhat  large  entropy  changes  around  the 
leading  edge,  where  the  H-grld  has  a  large  grid 
skewness  and  slope  discontinuity.1  These  numer¬ 
ical  errors  around  leading  edges  would  be  reduced 
If  a  C-grld  were  used  because  It  provides  a  better 
resolution  around  the  leading  edge.  A  standard 
C-grld  was  generated  for  a  core  turbine  stator 
vane.2  It  reveals  large  grid  skewing  on  the 
suction  side  of  the  flow  passage,  which  should 
Increase  numerical  errors.3*4  For  a  typical  tur¬ 
bine  blade,  no  single  grid  offers  satisfactory  grid 
properties  In  the  entire  turbine  stator  or  rotor 
passage. 

A  two-dimensional  0/H  patched  grid  used  In  a 
turbine  cascade  computation’  shows  the  slope 
discontinuity  at  the  0-H  grid  Interface  and  has 
extremely  large-aspect- ratio  meshes  near  the  lead¬ 
ing  edge.  The  numerical  scheme  used  In  that  work 
Is  a  cell-centered  scheme  based  on  the  Beam  and 
Warming  approximate  factorization.6  Although  the 
flexibility  of  the  scheme  Is  Indicated  In  two 
dimensions.  Its  likely  limitations  In  three  dimen¬ 
sions  suggest  an  alternative  approach. 

Two-  and  three-dimensional  composite  grids 
were  generated  In  an  attempt  to  Improve  the  grid 
quality  In  terms  of  mesh  control,  smoothness,  and 
orthogonality.  A  numerical  scheme  that  will  run 
on  this  composite  grid  Is  discussed. 

Composite  Grid 

A  stator  vane  ring  and  the  vane  geometry  at 
mean  section  are  shown  In  Figs.  1(a)  and  (b), 
respectively.  This  vane  ring  was  used  for  an  annu¬ 
lar  cascade  experiment  reported  In  Ref.  2.  To 
alleviate  numerical  errors  associated  with  the 
H-grld  skewness  for  the  blunt  leading  edge,  a 
standard  C-grld  (Fig.  2(a))  was  generated  to  exam¬ 
ine  any  possible  advantage  over  the  H-grld  for  this 
turbine  vane  geometry.  But  the  grid  became  very 
skewed  on  the  suction  side  of  the  passage  because 
the  blade  was  highly  cambered  (Fig.  2(b)).  In 
addition,  the  control  of  the  mesh  distribution  In 
the  blade  region  was  limited  with  the  C-grld, 
because  It  affected  the  mesh  distribution  In  the 
upstream  region.  Results  of  numerical  simulations 


1 


using  this  grid  are  likely  to  be  affected  by  the 
grid  skewing.3 


Numerical  Scheme 


Two-  and  three-dimensional  composite  grids 
were  generated  as  shown  In  Figs.  3(a)  and  (b), 
respectively.  The  blades  had  a  constant  profile 
from  hub  to  tip  and  were  stacked  at  the  trailing 
edge  (Fig.  3(c)).  Only  a  selected  number  of  grid 
lines  are  shown  for  Illustration.  This  composite 
grid  consists  of  the  C-grld  (or  0-grld)  In  the 
Immediate  vicinity  of  the  blades  and  the  H-grld  In 
the  upstream  region  and  In  the  middle  of  the  blade 
passage  between  the  C-grlds.  The  C-grld  (or 
0-grld)  portion  can  be  generated  by  using  either 
the  elliptic  method7  or  the  algebraic  method.3 
At  the  C-grld  and  H-grld  Interface  (or  0-H  Inter¬ 
face)  the  slope  continuity  was  preserved  so  that 
no  special  numerical  approximations  are  needed  for 
the  derivatives  at  the  Interface.  The  C-grld  (or 
0-grld)  Is  orthogonal  to  the  blade  surface  and  pro¬ 
vides  a  good  boundary  layer  resolution  near  the 
leading  edge.  This  composite  grid  has  better 
smoothness  and  orthogonality  and  provides  more 
flexibility  In  controlling  the  meshes  than  any 
single  grid  for  a  typical  turbine  blade  with  large 
camber  and  a  rounded  leading  edge. 

The  C-grld  portion  of  the  composite  grid  was 
generated  by  the  elliptic  grid  generation  code.3 
Two  subroutines  of  the  code  were  modified. 1 3  The 
outer  boundary  subroutine  was  modified,  to  Improve 
generation  of  the  periodic  boundaries,  by  using  the 
mean  camber  line  for  hlgh-solldlty  blade  rows.  The 
Inner  boundary  subroutine  was  modified  to  allow  a 
more  general  clustering  of  points  about  the  leading 
and  trailing  edges  of  highly  cambered  turbine 
blades.  Figure  4  Illustrates  the  portion  of  the 
C-grld  between  n  *  0  and  n  -  oc  that  was 
retained  for  the  composite  grid.  The  grid  lines 
at  n  =  0  and  n  *  nmax  are  Inner  and  outer 
boundaries,  respectively.  The  choice  of  nc  Is 
arbitrary  and  may  depend  on  the  extent  of  the  shear 
flow  region.  The  n'  coordinate  Is  a  translation 
of  the  n  coordinate  by  one  pitch  In  clrcumferen 
1 1  a  1  direction.  For  the  H-grld  In  the  middle 
between  the  blades,  the  C-grld  points  at  n  =  n 

and  n  =  n  ^  on  the  suction  side  and  at 
n1  =  and  n'  =  n|!fl  on  the  pressure  side  are 
used  In  the  cubic  spline  to  preserve  slope  contl 
nulty  at  the  C-H  interface  and  to  effect  a  smooth 
change  In  mesh  size.  The  three-dimensional  grid 
In  Fig.  3(b)  was  constructed  algebraically  from  the 
composite  grids  of  the  hub  and  shroud  surfaces  by 
using  the  cubic  spline  and  exponential  stretching. 

A  two  dimensional  composite  grid  was  also  gen 
erated  for  the  first- stage  stator  and  rotor  of  the 
space  shuttle  main  engine  (SSMt)  fuel  turbopump  tur 
bine,  figure  5  shows  the  construction  of  the  C  H 
composite  grid  for  the  rotor  blade.  Figure  6  shows 
the  composite  grid  generated  for  the  stator  vane. 
One  merit  of  the  composite  grid  Is  the  flexibility 
In  constructing  the  H  grid  In  the  upstream  region. 
One  can  choose  the  number  and  distribution  of 
meshes,  and  the  H  grid  can  be  extended  In  the 
upstream  direction  Independently  of  the  C  grid  part 
of  the  composite  grid  around  the  blade.  With  this 
flexibility  the  downstream  end  of  the  grid  for  the 


Semidiscrete  Finite-Volume  Scheme 

The  two-dimensional  C-H  composite  grid  has  an 
unusual  grid  point  that  Is  connected  to  more  than 
four  neighboring  nodes  (Fig.  8(a)).  At  this  point 
the  usual  differencing  techniques  cannot  be  applied. 
If  a  standard  finite-difference  scheme  were  used  to 
solve  a  flow  problem  on  this  composite  grid,  It 
would  be  difficult  to  treat  this  special  point.  A 
finite-volume  scheme1 1  to  be  used  with  this  com¬ 
posite  grid  presents  no  problem  and  requires  no 
special  treatment  because  each  Interior  cell  of 
this  composite  grid  has  only  four  neighboring  cells 
In  two  dimensions  and  six  neighboring  cells  In 
three  dimensions  (Fig.  8(b)).  The  f Inlte-volume 
scheme  Is  described  briefly  here. 

The  Euler  equations  In  Integral  form  can  be 
written  as 

it  $  "  dK  *  fff-  - 0  n) 

for  a  fixed  region  it  with  boundary  an.  Here 
w  represents  the  conserved  quantity,  F  Is  the 
corresponding  flux  term,  and  t  Is  time. 

A  convenient  way  to  ensure  a  steady  state 
solution  Independent  of  the  time  step  Is  to  sepa¬ 
rate  the  space  and  time  discretization  procedures. 

In  the  semidiscrete  finite-volume  scheme  one  begins 
by  applying  a  semidiscretization  In  which  only  the 
spatial  derivatives  are  approximated.  To  derive  a 
semidiscrete  model  that  can  be  used  to  treat  com¬ 
plex  geometric  domains,  the  computational  domain 
Is  divided  Into  quadrilateral  cells.  Assuming  that 
the  dependent  variables  are  known  at  the  center  of 
each  cell,  a  system  of  ordinary  differential  equa 
tlons  Is  obtained  by  applying  equation  (1)  sepa 
rately  to  each  cell.  These  equations  have  the  form 


dt  (S1jw1j)  *  Q1j 


where  j  Is  the  cell  area  and  Qu  Is  the 
net  flux  out  of  the  cell.  This  can  be  evaluated  as 
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£  u*k  fk  -  axk  v  <3> 

k-l 

where  f|(  and  g^  denote  values  of  the  flux 
vectors  f  and  g  on  the  kth  edge,  AX|<  and 
Ay^  are  the  Increments  of  x  and  y  along  the 
edge  with  appropriate  signs,  and  the  sum  Is  over 
the  four  sides  of  the  cell.  The  flux  vectors  are 
evaluated,  for  example,  by  averaging  the  values  In 
the  cells  on  either  side  of  the  edge: 


The  scheme  constructed  In  this  manner  reduces  to  a 
central  difference  scheme  on  a  Cartesian  grid  and 
Is  second  order  accurate  In  space  provided  that  the 
mesh  Is  smooth  enough.  It  also  has  the  property 
that  uniform  flow  Is  an  exact  solution  of  the  dif¬ 
ference  equations. 


stator  can  easily  be  matched  with  the  upstream  end  LU  Implicit  Scheme 
of  the  grid  for  the  rotor  In  the  vaneless  space 

between  the  stator  and  rotor.  Fiqure  7  Illustrates  The  composite  grid  Is  generated  to  be  used  for 

that  the  composite  grid  can  be  easily  constructed  both  Invlscld  and  viscous  flow  calculations.  For 
for  a  stator/rotor  Interaction  study. 


viscous  calculations  the  C-  or  0-meshes  must  be 
very  fine  to  resolve  the  boundary  layer.  And  It 
Is  likely  that  the  time  step  Imposed  by  an  explicit 
stability  bound  will  be  much  less  than  that  Imposed 
by  the  accuracy  bound  of  an  Implicit  scheme.  Since 
an  obvious  way  to  accelerate  convergence  to  a 
steady  state  Is  to  Increase  the  size  of  the  time 
step,  an  Implicit  scheme  Is  expected  to  have  a 
faster  convergence.  Although  the  alternating 
direction  Implicit  (A01)  scheme  has  been  valuable 
In  two-dimensional  problems.  Its  Inherent  limita¬ 
tions  In  three  dimensions  suggest  an  alternative 
approach.  An  LU  Implicit  scheme  was  demonstrated 
to  be  efficient  and  robust  for  both  external  and 
Internal  flows  In  a  broad  flow  regime. 1 2- 1 4  xhe 
scheme  was  extended  for  three-dimensional  flows  on 
an  H-grld.1^  This  scheme  will  be  used  for  the 
numerical  simulation  of  turbine  flows  on  the 
composite  grid. 


accuracy  drops  to  first  order.  The  unfactored 
Implicit  scheme  (eq.  (7))  produces  a  large  block- 
banded  matrix,  which  Is  very  costly  to  Invert  and 
requires  huge  storage.  An  unconditionally  stable 
Implicit  scheme  that  has  error  terms  at  most  of 
order  ( at ) ^  In  any  number  of  space  dimensions 
can  be  derived  by  LU  factorization 

ji  ♦  b  at  (dV  +  dV)J  ji  «■  s  At  (dV  ♦  oV)j 

x  AW  +  At  R  =  0  (8) 

where  Dx  and  Dy  are  backward  difference  oper¬ 
ators  and  DJ  and  DJ  are  forward  difference 
operators.  Ihe  reason  for  splitting  Is  to  ensure 
the  diagonal  dominance  of  lower  and  upper  factors 
as  well  as  to  make  use  of  the  built-in  Implicit 
dissipation. 


The  conservation  law  form  of  the  Euler  equa 
tlons  In  Cartesian  coordinates  for  two  dimensional 
flow  Is 


aw 

3F 

3G 

at 

*  ax 

♦  — 
ay 

<«> 


where  W  Is  the  vector  of  dependent  variables  and 
F  and  G  are  convective  flux  vectors: 


W  * 


(p,pu,pv,pE) 


Here  A*,  A  ,  B*,  and  B  are  constructed  so 
that  the  eigenvalues  of  matrices  are  nonnega 
tlve  and  those  of  matrices  are  nonpositive. 

A+  «  \  (A  ♦  r A 1 ) ,  A"  =  1  (A  -  r A I ) 

B*  =  \  (B  ♦  r B 1 ) ,  B*  =  \  (B  -  rBI) 

where 


F  =  [pu,pu‘  ♦  p,pvu,u(pE  ♦  p) 


(5) 


r A  >  max( |Xft| ), 


r B  >  max(  | | ) 


2  T 

G  =  [pv.puv.pv  <■  p,v(pt  ♦  p)  ] 

where  p,  u,  v,  E.  and  p  are  density,  velocity 
components,  total  energy,  and  pressure.  The  pres 
sure  Is  obtained  from  the  equation  of  state 

P  =  p(T  -  l)|t  -  \  (u2  ♦  v2)j  (6) 

where  y  Is  the  ratio  of  specific  heats.  These 
equations  are  to  be  solved  for  a  steady  state 
aw/at  =  0.  where  t  denotes  time. 

Let  the  Jacobian  matrices  be 

af  aG 


and  let  the  correction  be 


where  n  denotes  the  time  level. 

The  linearized  Implicit  scheme  for  a  system 
of  nonlinear  hyperbolic  equations  such  as  the  Euler 
equations  can  be  formulated  as 

1  V  6  At  ( Dx  A  ►  DyB)j  AM  *  At  R  -  0  (7) 

where  I  Is  the  Identity  matrix  and  R  Is  the 
residual 


Here  X/*  and  kg  represent  eigenvalues  of 
Jacobian  matrices.  Equation  (8)  can  be  Inverted 
In  two  steps.  The  LU  Implicit  scheme  needs  the 
Inversion  of  sparse  triangular  matrices,  which  can 
be  done  efficiently  without  using  large  storage. 
This  scheme  has  only  two  factors  In  three 
dimensions . 

C one  l.udlng  Rema  rk s 

Iwo  and  three-dimensional  composite  grids 
were  generated  to  Improve  grid  quality  for  an  annu 
lar  turbine  cascade  In  terms  of  smoothness,  reso¬ 
lution,  and  orthogonality.  Ibis  composite  grid, 
which  has  a  C  type  (or  0  type)  grid  In  the  Immedl 
ate  vicinity  of  the  turbine  blade,  provides  a  good 
boundary  layer  resolution  around  the  leading- edge 
region  for  viscous  calculation,  has  orthogonality 
at  the  blade  surface  and  slope  continuity  at  the 
C-H  (or  0-H)  Interface,  and  controls  the  mesh  dis¬ 
tribution  In  the  upstream  region  without  using 
excessive  grid  points.  This  composite  grid  elimi¬ 
nates  the  undesirable  qualities  of  a  single  grid 
when  generated  for  a  typical  turbine  geometry. 

The  C  H  composite  grid  has  flexibility  In  con 
structlng  the  H  grid  In  the  upstream  region  so  that 
the  grids  for  the  stator  vanes  and  rotor  blades  of 
the  SSME  fuel  turbopump  turbine  can  easily  be 
matched  In  the  vaneless  space  and  can  be  used  for 
the  stator  rotor  Interaction  study. 


R  =  DxF(Hn>  ♦  DyG(Wn) 

Here  0X  and  Dy  are  central  difference  oper 
ators  that  approximate  d/ax  and  a/ay. 


A  finite  volume  lower  upper  (LU)  Implicit 
scheme  Is  to  be  used  In  solving  the  turbine  flows 
on  the  composite  grid.  This  grid  has  a  special 
grid  node  that  Is  connected  to  more  than  four 
neighboring  nodes  In  two  dimensions  and  to  more 
than  six  nodes  In  three  dimensions.  But  the 
finite  vo'ume  approach  poses  no  problem  at  the 
special  r  ilnt  because  each  Interior  cell  has  only 
four  neighboring  cells  In  two  dimensions  and  only 


If  B  =  1/2,  the  scheme  remains  second  order 
accurate  In  time;  for  other  values  of  B  the  time 


six  cells  In  three  dimensions.  The  finite-volume 
LU  Implicit  scheme  was  proven  to  be  robust  and 
efficient  In  a  broad  flow  regime  and  Is  expected 
to  yield  accurate  solutions  on  the  Improved  compos 
1 te  grid. 
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